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THE FOURTH MOMENT OF DIRICHLET L-FUNCTIONS 



MATTHEW P. YOUNG 

Abstract. We compute the fourth moment of Duichlet L-functions at the central point 
for prime moduh, with a power savings in the error term. 



o 

^ in large part to extensive applications. Yet, these moments are seen to be natural objects 

^ to study in their own right as they illuminate structure of the family and display beautiful 

* 1 /" 



1. Introduction 

Estimating moments of families of L-functions is a central problem in number theory due 



symmetries. 

The Riemann zeta function has by far garnered the most attention from researchers. 
Ingham [In] proved the asymptotic formula 



^ I led + ^t)\'dt = a,{logT)' + OiilogTf 



c5. T 

where 04 = (27r^)~^. Heath-Brown |H-B2j improved this result by obtaining 

^ (1.1) U |C(| + zt)rrft = X^a,(logT)^ + 0(T-i+^) 

in i=o 

^ . for certain explicitly computable constants (see (5.1.4) of |CFKRS] ). Obtaining a power 
O i savings in the error term was a significant challenge because it requires a difficult analysis 
^ ' of off-diagonal terms which contribute to the lower-order terms in the asymptotic formula. 
The fourth moment problem is related to the problem of estimating 

:a 



(1.2) Y,d{n)d{n + f) 

n<x 

uniformly for / as large as possible with respect to x. Extensive discussion of this binary 
divisor problem can be found in |M1] . The sum (11. 2p can be transformed into a problem 
^ ■ involving Kloosterman sums. The strength of Heath-Brown's result (11. ip relies on the Weil 
bound. Using the spectral theory of automorphic forms (the Kuznetsov formula), Iwaniec 
( |Iwl] . Theorem 3) showed 

(1.3) / |C(i + ^t)|"c/t<T2/3+^ 



Notice that Heath-Brown's result gives fOjl with replaced by T^/^ With some extra 
work, these mean-value estimates can be turned into subconvexity estimates for the zeta 
function (although Weyl's method already gives C(| + it) < with far less effort). 

Motohashi has proved a beautiful exact formula for a smoothed version of the fourth 
moment of zeta in terms of the cubes of the central values of all automorphic (degree 2, level 
1) L-functions ( [M2] , Theorem 4.2). The result is difficult to describe precisely because the 
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dependence on the smoothing function involves elaborate integral transforms. Motohashi 
considers 

/oo 
C(i + a + it)C{l + (3 + tt)al + 7 - it)al + 5 - tt)g{t/T)dt, 
-oo 

where initially the parameters a,/3,7,5 have large real parts, and uses the Kuznetsov for- 
mula to develop a meromorphic continuation of Mg(a, /3, 7, 5) to include the origin. These 
parameters (or shifts) are helpful for developing the main term in the asymptotic formula 
for the fourth moment because certain residue computations become simplified (e.g., it is 
much easier to compute the residues of five simple poles than it is to compute a residue at a 
fifth order pole). The smoothing function g can be replaced with the sharp truncation t <T 
using standard techniques (however, unsurprisingly, the error term becomes larger). 

Recently, five authors |CFKRS] have developed a recipe for conjecturing moments of fam- 
ilies of L-functions, including all lower order terms. The presence of these shifts plays an 
important role in their conjectures. Besides allowing for simpler residue computations, the 
presence of these parameters makes symmetries of the family visible. For example, Mg is 
invariant under switching a and /3 or 7 and 6. Application of the functional equation for 
zeta shows relations under a — )■ —7, 7 — t- —a, and other similar relations. The shifts are 
useful too since for example one can differentiate with respect to them to study moments of 
the derivatives of a family. 

The family of all primitive Dirichlet L-functions of modulus q is analogous in some ways 
to the Riemann zeta function in t- aspect. However, there are significant differences which 
cause the family of Dirichlet L-functions to be more difficult to study. Ramachandra [R] 
conjectured that for prime p, 

p-' \Lihx)\'-i2n'r\^ogp)\ 

x{mod p) 

Heath-Brown |H-B1] proved 

(I'S) ^ r i^(i-^)i'-^n^^a°^^)'+o(2"<..^„og,)-'). 

x{modg) p\q 

where 4'*{q) is the number of primitive characters modulo q, u}{q) is the number of distinct 
prime factors of q, and the sum is over all primitive characters modulo q. For almost all 
q this is an asymptotic formula, however if uj{q) > (log log g) / log 2 the error term is not 
smaller than the main term. Recently, Soundararajan |Slj improved the error term in fll.Sp 
so that an asymptotic formula does indeed hold for all q. 

In this paper we obtain the asymptotic formula with a power savings for prime moduli. 

Theorem 1.1. For prime q ^ 2, we have 

(1-6) \L{hx)\' = i2^^{logqy + 0{q-^-^n, 

for certain computable absolute constants Ci. Here —5/512 = —1/80 + ^/40, where 9 = 7/64 
is the current best-known bound on the size of the Hecke eigenvalue X{n) of a Maass form, 
that is, a bound of the form \X{n)\ < d{n)n^ . 

This result is the analog of f lLip for Dirichlet L-functions of prime moduli. The method 
of proof can possibly generalize to handle certain variations on fll.6p with twists of fixed 
Dirichlet L-functions. Such variations could be used to prove simultaneous nonvanishing 
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results for four Dirichlet L-functions. For this type of variation it is imperative to have a 
power saving in the error term because the main term may be as small as a constant. 

The study of ([LB]) with L(l/2, x)^ replaced by L{f x x, 1/2) for / a Hecke cusp form has 
similarities to the analysis in this paper. The difference amounts to replacing the divisor 
function with the Hecke eigenvalues of the cusp form. In this paper we crucially use the fact 
that the divisor function can be written as 1 t^t 1. It is a difficult and interesting problem to 
obtain a power savings in the second moment of Dirichlet twists of a Hecke cusp form. 

In order to obtain power savings in (II. 6p it is necessary to obtain an asymptotic formula for 
the off-diagonal terms. The results of Heath-Brown and Soundararjan arise from bounding 
the contribution of these off-diagonal terms. To elaborate, the problem of estimating (11.61) 
essentially reduces to the analysis of the following divisor sum 

d\q m=n (mod d) 

{mn,q)=l 

where V{x) is a smooth function with rapid decay, satisfying V{x) ~ 1 for x small. The 
contribution of the diagonal terms m = n in the sum is easily computed and gives the main 
term of (11.51) . 

The off-diagonal terms m ^ n are much more difficult to analyze. One of the primary 
difficulties in treating (II. 7p is the large ranges of summation of m and n. Consider the 
sum Bm,n consisting of those summands in (11.71) in the dyadic segments M < m < 2M, 
N < n < 2N, where MN < (by symmetry suppose > M). If M x 1 then N can be 
as large as but if M x g then N ^ q also. It is not surprising that different techniques 
are required to estimate Bm,n for these different regions. In fact, the relative sizes of M and 

has a large effect on how large Bm,n is. Qualitatively, Bm,n is small if M and are far 
from each other, but when M and N are close to the same size, then Bm,n contributes to 
the main term of (11.61) . To see this former assertion heuristically, we use the fact that the 
divisor function is uniformly distributed in arithmetic progressions, that is if (m, q) = 1 

(1.8) ^^'') = 0rT + (error). 

n<x n<x 
n=m (mod q) {n,q)=l 

It is known unconditionally that (11.81) holds for q < xi"^ with an error of size ©(x^"*"^), 
following from Weil's bound for Kloosterman sums; see Corollary 1 of |H-B2j for a proof. 
Note that if q is prime, then 

V^"^^^ M<m<2M N<n<2N N<n<2N 

{m,q)=l n=m (mod q) {n,q)=l 

which is small provided the error term of (11.81) is smaller than the main term. 

Improving the uniformity for which (II. 8p is true is a challenging open problem. It seems 
natural to expect that the divisor function is evenly distributed across arithmetic progressions 
for q < x^~^. Fouvry [FJ proved this is true on average over x"^^^^^ < q < x^~^ , and Fouvry 
and Iwaniec |FoIj have produced results which cover x'^^'^^^ < q < x'^^'^^^ for special values 
of q . These results fix the arithmetic progression m (mod q) and average over g; for the 
application of estimating (II. 7p . q is fixed but m is allowed to vary. One of the main difficulties 

1 3 

in this work is treating the range of summation in (I1.7P where M x ga, x gz, which is at 
the edge of the range where the error term of 0(x3^^) for (II. 8p is barely insufficient. 
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The expectation fll.Sp breaks down when x and q are of comparable size. Instead, main 
terms are formed from the summands with M ^ N ^ q. To see this heuristically, note that 
in this range the sum (11 .7^ is mimicked by 



which can be seen by opening the square and summing over a. Clearly this expression cannot 
be small because the sum of m ^ q, m = a (mod q) essentially picks up one term which 
certainly cannot approximate the average behavior of the divisor function! 

This dichotomy is somewhat analogous to the (smoothed) fourth moment of the zeta func- 
tion, where averaging over t aspect forces m and n to be close to each other (because there 
is an integral of the form J{n/my^g{t)dt, which is small unless m and n are close). It is 
remarkable that the same phenomenon occurs for Dirichlet L-functions (the main contribu- 
tion coming from m and n close) but for completely different reasons. If we combined our 
average over characters with a short t-aspect integration then our problem would become 
much simpler; however, this would destroy any nonvanishing applications. 

One apparent difficulty in finding an asymptotic formula for (11. 7p is that the values of 
V{mn/q^) when m and n are around q are at the transition range of V . However, this type 
of behavior has been encountered by Soundararajan ([22], discussion following (5.16)) and 
Kowalski ([K], p. 155) in their studies of other families of L-functions. The analysis of the 
main term in this region of summation leads to a contour integral which can be computed 
exactly using symmetry properties of the integrand, which in turn relies upon the functional 
equation for the Riemann zeta function. In our work this occurs in the proof of Lemma 18.31 

In this paper we use essentially two different methods of handling the sum (11.71) . depending 
on the sizes of M and A^. When is significantly larger than M, then we treat the sum 
over n as the divisor function in arithmetic progression. Using only Weil's bound in a 
straightforward way, we could obtain the necessary asymptotic formula for N/M > 
We succeeded at extending the range of summation of n to smaller values by exploiting 
extra savings by averaging over m. These arguments are presented in Section HI 

The region where M and N are relatively close is treated using similar techinques as 
Motohashi in his work on the fourth moment of the Riemann zeta function |M2j . The 
Kuznetsov formula plays the key role in estimating the binary divisor sum (II. 2p : the range of 
uniformity of / with respect to x depends on the best-known bound on the size of the Hecke 
eigenvalues of Maass forms (see |Mlj . Theorem 5). Assuming the Ramanujan-Petersson 
conjecture would provide an asymptotic formula for (ll.7p in the region N/M < 

Actually, our treatment is different than Motohashi's because we chose to use approximate 
functional equations in order to have formulas a priori valid in our region of interest (e.g., 
at the central point of the critical strip). This feature causes the main terms to be captured 
in a more straightforward way in our work. In addition, approximate functional equations 
explicitly display how large variables are with respect to each other. Nevertheless, the basic 
strategy closely follows Motohashi. 

One curious aspect of our work is an approximate generalization of Motohashi's formula 
for the fourth moment of the zeta function in terms of the third moment of Hecke-Maass 
L-functions. We arrive at a certain formula involving products of three Hecke-Maass L- 
functions twisted by the g-th Hecke eigenvalue; see (19.121) below. An important difference 
between these formulas is that our (19.121) is only helpful when M and A^ are restricted to be 
fairly close. 



a (mod q) 




moiq 
m=a (mod q) 
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We have chosen to compute the shifted fourth moment of Dirichlet L-functions, that is 
we include the parameters a, /3, 7, 5 similarly to f ll.4p . Doing so allows for a more structural 
viewpoint of the main terms, and also allows us to check the conjecture of the five authors. 

1.1. Acknowledgements. I am most grateful to Brian Conrey for extensive discussions on 
this material and for suggesting this problem. I benefitted from conversations with David 
Farmer, Peng Gao, and Soundararajan. I especially want to thank Henryk Iwaniec for a 
number of useful discussions. 



2. Notation and background 

2.1. Dirichlet L-functions. In this section we briefly recall some standard facts needed 
about Dirichlet L-functions. Let q be an odd prime (g > 3), x be a primitive character 
modulo g, and be the Dirichlet L-function L(s,x) = Xln ^('^) "^^^ completed 

L-function satisfies the functional equation 

A(., x) = ) ' r x) = e(x)A(l - s,x), 

where r(x) = Y.x (mod <?) x{x)e j is the Gauss sum, and 




In its asymmetric form, the functional equation reads 

L{s,x)=e{x)X{s)L{l-s,x), 

where 



1 

-1. 



(2.1) X{l + u)=(^ 



-„ r 



g\ " \ 2 



2 — 11+n 



71/ Y I 2+"+° 



The central quantity of interest in this paper is 

(2.2) M(a, /5, 7, S) = ^ L{\ + a, x)L{\ + x)L{\ + 7, x), L{\ + <5, x), 

'x (mod q) 

where the + indicates the summation is restricted to primitive even characters and 0*(g) is 
the number of primitive (odd or even) characters. Throughout this work we assume that the 
shifts a,/3,7,(5 are all sufficiently small with respect to q (say, <C (logg)~^). It is natural to 
split the family separately into even characters and odd characters because the two families 
have different gamma factors in their functional equations. In this work we concentrate 
almost exclusively on the even characters because the case of the odd characters is similar 
(we could treat both cases simultaneously but it would clutter the notation). In Section [5751 
we briefly describe the necessary changes to treat the odd characters. 
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2.2. The conjecture for the fourth moment. 

Conjecture 2.1. [CFKRSj For any q ^2 (mod 4), and with shifts <^ (logg)~^, 

C,(l + a + 7)C,(1 + a + 5)C,(1 + /3 + 7)C,(1 + /? + 5) 



(2.3) M(a,/3,7,5) 



Cg(2 + « + /? + 7 + 5) 

C,(l - « + /3)Cg(l - g - 7)C,(1 + /9 + '^)Cg(l - 7 + ^) 
C,(2-a + /3-7 + 5) 

C,(l - g + /3)C,(1 -a- 6Xg{l + /3 + 7)Cg(l + 7 - '^) 
C,(2-a + /3 + 7-5) 

C,(l + a - /3)Cg(l + + ^)C,(1 -P- 7)C,(1 - 7 + '^) 
C,(2 + «-/3-7 + 5) 

Cg(l + a- /3)C,(1 + + 7)C,(1 - /g - '^)C,(1 + 7 - ^) 
C,(2 + a - /3 + 7 - 5) 



^ y C,(l - g - 7)C.(1 - c:^ - '^)C,(1 - /3 - 7)C.(1 - /3 - (^) ^ -1/2+.^ 
+ ^"''^'^'^ C,(2-a-/3-7-5) + ^' 

w/iere C<7(s) = C(s) Ilpigll - P^'') o-i^d Xa,y = X(i + a)X{^ + 7), eic. 

In the case of q prime one may replace each occurence of (q with ( without ahering the 
error term. There are no primitive characters modulo g if g = 2 (mod 4). 

The main terms on the right hand side clearly exhibit the same symmetries as the mo- 
ment must (arising from trivial permutations of the variables as well as applications of the 
functional equation of the Dirichlet L-functions). A similar conjecture holds for the odd 
characters; the only change is to have the X factors given by ( 12. ip depend on the parity. 

It is important not to underestimate the psychological advantage of having the main 
terms produced ahead of time. The actual computation of these main terms is by no means 
straightforward; there are many cancellations and combinations of various terms that occur 
throughout the work. Although it may appear that (12.31) is a messy expression, in fact it 
arises after many simplifications. 



Theorem 2.2. For q > 3 prime, Coniecture \2. 1\ holds hut with an error of size 0{q 



-5/512H 



Theorem 11.11 follows by taking the limit as the shifts go to zero. It is not obvious at first 
glance that the main term on the right hand side of (12. 3p is holomorphic with respect to the 
shift parameters, but this is proven in a more general setting in |CFKRS] : a computer can 
also easily verify it in this case. Of course the left hand side is holomorphic in the shifts. 

In the course of the work we assume that each of the shifts lies in a fixed annulus with 
inner and outer radii x (logg)~^; we further suppose that these annuli are separated enough 
so that any linear combination of the form eia + e2/3 + €37 + €46 with each ej G {—1, 0, 1}, 
not all zero, is 3> (logg)~^. We initially prove Theorem 12.21 with this restriction in place; by 
the holomorphy of the left hand side of (12.31) . as well as the combination of main terms on 
the right hand side of (12. 3p . we see that the error term also must be holomorphic in terms 
of the shifts. By the maximum modulus principle, we can then extend the result to all shifts 
<^ (logg)^^ with the same error term. 



2.3. The orthogonality formula. 
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Proposition 2.3. // {ab, q) = I then 

(2.4) E^^(«M^) = ^ E <f>(d)Md). 

X (mod g) d\{q,a±b) 

The sum on the left hand side vanishes if {ab, q) ^ 1. The condition d\a±b should be taken 
with multiplicity (i.e., if d\{a + b,a — b), it is counted twice). 

The proof is standard and appears in many sources such as |H-Blj . |Slj . The odd charac- 
ters have a similar orthogonality relation except the terms with d\a + b are subtracted from 
the terms with d\a — b. Also, note that for q prime, = q — 2, and the number of even, 

primitive characters is |(g — 3). 

2.4. Approximate functional equation. To get a useful formula for M{a,(3,'y,6) we 

use an approximate functional equation. There are a variety of choices in how to represent 
|L(|, we have chosen to use the functional equation of L(s, x)^L{s, x)^ to write \L{j, 
as the sum of two sums of length approximately g^. Actually, we need a formula for + 
a,x) ■ ■ ■ + S,x), which prevents us from easily lifting a formula from, say. Theorem 5.3 
of IIKI. however the derivation of the formula follows standard lines. We have 



Proposition 2.4 (Approximate functional equation). Let G{s) be an even, entire function 
of exponential decay as \s\ — )■ oo m any fixed strip \Re{s)\ < C and let 

(2-5) Vo,,i3,^,s{x) = — I —ga,i3,^A^)^~'ds, 

where 

(2.6) 9a,t3,^As) = TT-'^ 



27rz ./(i) s 



Furthermore, set 

(2.7) X^,p,,,s = X(| + «)X(i + /3)X(i + 7)X(i + 6). 
Then 

(2.8) L(i + a, x)Lil + (3, x)L{\ + 7, x)L{\ + 5, x) 



E x(mim2)x(^i^2) ( ^i"^2»i^2 ^ 

i+« 1+^ 1+^ 1+5^/3.7,51 2 1 



\ - X("^i"^2)x('^i^2) ( mim2nin2\ 

+ -^a,/3,7,<5 2^ i_„ i_« i_5 -/3 -7 -5 ( 



2 ^ 2 ^ 2 ' „ 2 \ q / 



^2 ""1 ""2 

We impose additional conditions on G which we separate here for ease of reference. 

Definition 2.5. Let Qa,i3,'y,s{s) be an even polynomial in s with the following properties: it 
takes the value 1 at s = 0; it is rational in the shifts a,(3,'y,6; it is symmetric in the shifts; 
it is invariant under a — )■ —a, (3 — >■ — /3, etc.; it is zero at 2s = —a — 7, s = ^±a (as well 
as (3 + 6, I ± /9, etc., by symmetry) . Then set G{s) = Qa,i3,-y,s{s) exp(s^). 
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Proof of Proposition \2.J\ Let 

K^fiaA^) = A(| + s + a, x)A(| + s + /3, x)A(| + s + 7, x)A(i + s + 5, x) 
and consider 

h= Aa,/3,7,5(s) ds. 

J{1) s 

Move the line of integration to (—1), passing the pole at s = 0. Let I2 be the new integral. 
The residue at s = is 

Aa,/3,^,5(0) = A(i + a, x)A(i + 13, x)A(i + 7, x)A(i + 6,x)- 

After the change of variables s — )■ — s and the application of the functional equation 

Aa,/3,7,5( — s) = A_^^_5_Q, 

we obtain 

h = - A-^,^s,-a-p{s) ds. 

J(l) s 

Set 

^a,/3,7,<5(s) = L(i + a + s, + /3 + s, x)L{\ + 7 + 5, + ^ + -s, x), 



and let 

Then we have 



-t^a,/3,7,<5lUj - 7r~ ^a,l3,-f,5[S)— 7-7- dS 



2Ti 7(1) '• • • ' A4,j,-,,j(0) s 



An easy computation shows 



Expanding L^, „,(s) into absolutely convergent Dirichlet series and reversing the order of 
summation and integration completes the proof. □ 

2.5. Automorphic forms. We briefly summarize the material we require on automorphic 
forms in order to apply the Kuznetsov formula. See |M2] . |Iw3] . or jlK] for further details. 

Let Uj{z) be an orthonormal system of Maass cusp forms on S'L2(Z)\EI with Laplace 
eigenvalues j + k^. Each cusp form Uj{z) has the Fourier expansion 



= 2/^ ^ Pj{n)Ki^,i2'iT\n\y)e{ 



We assume that each Uj is an eigenfunction of the Hecke operators T„ defined by 



ad=n b (mod d) 
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with eigenvalues Xj{n). That is, T^Uj^z) = Xj{n)uj{z). By consideration of the reflection 
operator, we may furthermore assume that pj{—n) = ejPj{n). The Hecke eigenvalues are 
multiplicative and satisfy the Hecke relation 

(2.9) A,(m)A,H = J2 ^.•(^)' 

d\{m,n) 

In terms of Fourier coefficients, the Hecke relations give Pj{n) = pj{l)Xj{n). Estimating the 
size of an individual Hecke eigenvalue is an important problem. So far the best result is 

(2.10) |Aj(n)| < d(n)n^ 

with 9 = 7/64, due to Kim and Sarnak |KSj . The Ramanujan-Petersson conjecture asserts 
that 6* = is allowable. 

The Hecke-Maass L-function is initially defined by Lj{s) = Xj{n)n~'^, which converges 
absolutely for Re(s) > 1 by properties of the Rankin- Selberg convolution. Then Lj{s) 
continues to an entire function and satisfies the functional equation 

AM := .-r (t±3±i^^ r (i±ipa) i,,,) = .^A,(i - 

where Sj = (1 — ej)/2. The Hecke relations (12.91) translate to Lj{s) having the Euler product 

The Eisenstein series have the Fourier expansion 

(2.11) 7r-T(s)C(2s)E(z, s) = n-ris)a2s)y' + 7r^-^r(l - s)C(2 - 2s)y'-' 

+ 2y^ J2 \n\''^'(^i-2s{\n\)K^_i{27i\n\y)e{nx), 

where 

(2.12) ax{n) = J2d^- 

d\n 

The Fourier expansion furnishes the meromorphic continuation of E{z, s) to s G C. 

Let Uj^k, 1 < J < dim S'fc(S'L2(Z)) be a complete orthonormal Hecke basis of the classical 
weight k cusp forms. These have the Fourier expansion 

The Ramanujan-Petersson conjecture is known for holomorphic forms due to Deligne jP] . 

2.6. Kloosterman sums and Kuznetsov formula. The Kuznetsov formula is a trace 
formula relating sums of Kloosterman sums to Fourier coefficients of automorphic forms. 
This technology can show that there is considerable cancellation in the sum of Kloosterman 
sums. Furthermore, it provides a separation of variables of m and n in Y2c S{m, n; c) which 
is conducive to obtaining additional savings in summations over m and n; see [Iw2] and |DIj . 



10 MATTHEW P. YOUNG 

Theorem 2.6 (Kuznetsov formula). Let g be aC^ function satisfying g{0) = and g^^\x) <^ 
{x + l)"^""^, j = 0, 1, 2 and suppose m,n > 0. Then 

(2.13) c ^ \~^~) ^ 



fc=0 (mod 2) i 

Here Aig and Mg are the following integral transforms 



■^fl(^) = -NA- / Jk-i{x)g{x) 



[At^iY Jo X 
For the opposite sign case, 

(2.14) 2^ ' ' ' ^7 ( 1 =}_^Pj{m)pj{~n)Kg{Kj) 

c 3 



/C,(i) = 2 / A2.,(i)g(i)— . 



where 



For proofs we refer to Theorems 16.5 and 16.6 of |IKj or Theorems 2.3 and 2.4 of [M2] (we 
borrowed some notation from both sources). 

2.7. The Estermann D-function. The Estermann D-function is defined by 

(2.15) Dis,X,j) = J2^^hj)^ 

where {h, I) = 1. The analytic properties of D are useful for understanding the behavior of 
the divisor function. We have 

Lemma 2.7. For any fixed A G C, D{s, A, y) is meromorphic as a function of s, and satisfies 
the functional equation 

(2.16) D{\ + s, A, y) = 2(27r)-i-^+2.p( 1 _ ^^^^i ^ ^ _ ^)^a-2s 

X D(| - s, -A, y) cos(f ) + D(i - s, -A, -y) sin(7r(s - I)) 
// A 7^ ^/ien /ias simple simple poles at s = 1 and s = 1 + A with respective residues 

(2.17) ri+^C(i-A), ri-^c(i + A). 

We refer to |M2] . Lemma 3.7 for proofs. The functional equation of the Estermann function 
is essentially equivalent to the Voronoi summation formula (Theorem 4.10 of |IK] ) . 
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2.8. Conventions. We use the common convention in analytic number theory that e denotes 
an arbitrarily small positive quantity that may vary from line to line. Furthermore, a, /3, 7, 
and 5 are complex numbers that are sufficiently small in comparison to e (so for example we 
may say C(l + a + s) is holomorphic for Re(s) > e). 

3. The structure of the fourth moment 

3.1. Averaging the approximate functional equation. Using Propositions 12.31 and 12.41 

we may average the approximate functional equation to obtain a formula for the fourth 
moment M(a,/3,7,(5) defined by (12. 2p . Write 

(3.1) M(a, /3, 7, 5) = Ai,,(a, /3, 7, 6) + A_i,g(a, (3, 7, 6), 

where Ai^q is the contribution from the 'ffist part' of the approximate functional equation, 
and A_i^q is the second part. It suffices to compute Ai^g since 

(3.2) ^-l,g(a, 13, 7, = ^a,/3,7,<5^1,q(-"> "Z^, "7, S)- 

We have 

•4., - ^ E E E E E TTTHbTTTi^--' (^^^) ■ 

^ d\q mim2=±nin2 (mod d) m| Uf \ ^ / 

{mim2nin2 ,q)=l 

By changing variables m = mim2 and n = nin2, we may rewrite this as 



(3.3) EE^^M^w(^) 

d\g m=±n (mod d) 

{mn,q)=l 

Let Ai be the same sum as g but with the condition {mn, q) = 1 omitted. Since q is 
prime, a trivial estimation gives 

(3.4) A,, = Ai + 0(g-i+^). 

3.2. The diagonal terms. In this section we compute the diagonal contribution Ad of the 
terms m = n in Ai given by (13.31) . To be precise we are computing the diagonal contribution 
of the terms where m = n (mod d) and not those with m = —n (mod d). We compute 



^D = ^)2^^(d)Md) ^ ^^^^^ v^,p,,A-\ 

^ d\q (n,q) = l ^ 



(n,g)=l 

Now the problem reduces to a standard exercise in analytic number theory. We use the 
integral representation of V , that is (12. 5p . to write the sum as an integral of a Dirichlet series 
and develop the asymptotics by moving the line of integration. Precisely, we have 



(n,g)=l 



From the Ramanujan identity 

>p ax{n)a^{n) _ C(f )C(^ - A)C(f - i^)C,{v - \ - y) 
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we obtain 



Ad = — fl'a,/3,7,5(s)g^'^a,/3,7,5(s)cis, 

2m Jn) s 



where 



'(1) 

C(l + a + 7 + 2s)C(l + a + 6 + 2s)C(l + /9 + 7 + 2s)C(l + (3 + 6 + 2s) 



C(2 + a + /3 + 7 + 5 + 4s) 

Z has simple poles at 2s = —a — 7, —a — 6, — /3 — 7, and — /3 — 5 and is otherwise holomorphic 
for Re s > — 1/4 + e. Note that G{s) vanishes at these poles of Z by our choice in Definition 
12.51 We move the line of integration to the line —1/4 + e, passing the pole at s = only. 
We obtain = (Residue) + /, where I is the contribution from the new line of integration. 
We easily obtain 

|/| < g"5+^. 

The residue at s = gives 

3.6 Yi{a, (3; 7, 6) := — — — — — , 

C(2 + a + /3 + 7 + a) 

which is precisely the first term in Conjecture 12.11 (up to 0(g^^^^)). We summarize this 
calculation with the following 

Lemma 3.1. With G as in Definition \2.5l we have 

(3.7) Az, = Fi(«,/3;7,5) + 0(g-i/2+-), 
and similarly the contribution of the diagonal terms to A^\ is 

(3.8) A^u = Y.,{a, (3; 7, 6) + 0{q-"^^'). 
Here 

(3.9) y_i(a, P- 7, 5) = X,,;3,7,5>^i(-", -/S; -7, -8)- 

If we did not choose G to vanish at the poles of Z then there would be four more extra 
terms in each of (13.71) and (13. 8p . However, these extra terms would all cancel out in the final 
assembly of M(a, /3, 7, 5); this satisfying (but complicated) calculation is carried out in the 
first version of this article on the arxiv. 

3.3. Decomposition of the off-diagonal terms. Here we investigate the contribution 
Ao to Ai from the off-diagonal terms m = n (mod d), m ^ n. The treatment of the 'dual' 
terms Aq corresponding to m = —n (mod d) proceeds in much the same way but must be 
executed separately, and is carried out in Section [71 
We need to compute 

(-0) ^o^^^^m^il) EE;^^^^i:^w(^), 

a\g m=n (mod a) 

where the prime indicates that the case m = n is excluded. 

We break the sum Aq into two pieces depending on whether m < n or m > n and write 
Ao = Bo + B'q accordingly. By symmetry, it suffices to compute Bo = Bo{a, /3, 7, S), since 
B'Q^a, /3, 7, 6) = Boi'j, S, a, /3). We record the decomposition, 

(3.11) A,=Ad + Bo + B'o + Ao. 

where each quantity has parameters a, /3, 7, 6 (in that order). 
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We require a partition of unity with some special properties. To begin, we take a partition 
of unity {Wm{x)} on M+ where each Wm{x) has support in the dyadic interval [M, 2M], and 
furthermore Wm{x) = W{jj) where W is a fixed smooth, compactly supported function. 
Here M varies over a set of positive real numbers, with the number of such M less than 
X being O(logX). Such a partition can be constructed very easily. We then create the 
partition of unity {WM,N{x,y)} on IR+ x IR+ by taking products of the Wm{x). 

We apply this partition of unity to Bp and write Bp = J2m n Bm,n, where Bm,n is the 
same expression as Bq but weighted with WM,N{^,n). Since m < n we may assume 

(3.12) M <^N, 

a convention that holds throughout this paper. Due to the rapid decay of V{x), we may 
assume MN < 



3.4. Overview of the proof. We write 

(3.13) Bm,n = (Main term)M,Ar + Em,n, 

for a certain main term that we do not explicitly write here due to its complexity. The 
discussion of the main terms takes place primarily in Sections |6] and [HI In Section 16.21 we 
prove that the size of the main term as a function of M and is given by 



(3.14) {Main term) M,N <t: M 2 N- 2 q' 



so it is a bit of a misnomer to call this a 'main term' when M and are not near each other. 

In Section H] we prove the following estimate for Bm,n that is applicable when N is signif- 
icantly larger than M. 

Theorem 3.2. ForM-^N, MN < 5^+^ we have 

(3.15) Em,n < (M^N-^ + q-iN'-i + g^A^"^M~A)g^ 

Actually we prove the bound fl3.15p but with Bm^n replacing Em,n- Notice the first term 
above is larger than the bound fl3.14p for the 'main term', so in fact the two bounds are 
equivalent. We do not attempt to extract a main term in this analysis, but still safely claim 
the same main term relevant in the range where M and A^ are close. 

To cover the range where M and A^ are fairly close, we prove the following in Section [Hi 

Theorem 3.3. With 9 = 7/64, M <^ N , and MN < we have 

(3.16) Em,7v < g"^+^+^M-5A^i 

The various main terms combine in a rather complicated way. When summed over all M 
and A^, all the main terms from Aq, Aq, etc. add up to form the quantity in Conjecture 12.11 
For some small values of MN we can do no better than the trivial bound 

Vmn , 

(3.17) Bm,n < q'. 



Note again that this bound on Bm,n is larger than the bound f l3.14p for the main term 
provided N > q (which is trivially satisfied since Bm,n is void otherwise), so this translates 
into a bound on Em,n- 
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Deducing Theorem 12.21 from the above results is then an exercise in finding the maximum 
of a piecewise linear function. Precisely, let 

Li(a, 6) = -1 + ifl + L2{a, b) = + 9 - + ^b, 

11 11 3 2 1 

L^fa, b) = -a 6, LJa, b) = b, L^(a, b) = a b. 

■i\ ^ ) 2 4 ' 4^ ' ^ 4 4 ' ' ^ 10 15 6 



We require 



max max min(Li(a, 6), -^2(0, &), I/j(a, b)). 

ie{3,4,5} 0<a<6 
a+b<2 



Along the line a = the maximum value is — ^ at 6 = ||. When a = b the maximum is 
— I + 6' and along the line a + b = 2 the maximum is + ^ (although if 9 was not known to 
be less than 6/43 then — 3^ + y would be the maximum here). The maximum at an interior 
point must occur when Li = L2 = Li for i = 3, 4, or 5. The case i = 3 is a=^ + 9,b = ^, 
which gives the bound —-^ + |. The case z = 4 gives — ^ + | from a = ^ + 9,b=^ — ^9. 
The case i = 5 has a = ^ + 9,b = ^ — ^9, and gives —-^ + {q, completing the proof. Note 
that our method requires the bound 9 < 1/6 which is very deep. 

4. The divisor function in arithmetic progressions 

In this section we present a handful of different estimates that taken together prove The- 
orem 13.21 We largely use classical techniques of analytic number theory such as Poisson 
summation, Cauchy's inequality, 'gluing' of variables to create a longer variable, reciprocity 
laws, the Weil bound for Kloosterman sums, estimates for exponential sums, etc. 

Solely for notational convenience we shall work with the case a = /3 = 7 = 5 = 0; the 
arguments extend easily to nonzero parameter values. 

4.1. Initial cleaning and statement of results. Recall 

d\g »n=n (mod d) 

Throughout this section, W{x) will denote a smooth function with support in a dyadic 
interval, which may vary from line to line (simply to avoid cluttering the notation). Since 
we treat the range M x in a different way, we assume N > cM for a sufficiently large 
constant c so that the term m = n is avoided by the support of the weight functions in (14. ip . 
As a simple first step, we separate the variables m and n in the test function V by using its 
Mellin transform. We have 

I f f q' VGis) 



2-ni ./(e) V MN ) s 



' Y^m^i) EE 

~~ ' m2n2 \ivi/ \i\ / 

a\q m=n (mod a) 



where Ws{x) = x ^W{x). Note that the dependence of Ws on s is very mild: 



x'—Wsix) < a;-^^(^)|R(s)| max 



x'—Wix) 
dx' ^ ' 



where Pj{s) is a degree j polynomial in s. In effect, this separation of variables has no cost 
since losses of size {MNy are absorbed by the factor q'^ in the bound (13.151) . and because 
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s ^G{s)g{s) has exponential decay in the imaginary direction. Thus, the problem reduces 
to bounding -B^/tv which is defined to be the same sum as Bj^i ^ but with V removed. Let 

^M,N{'m) = 



d\q n=m (mod ci) 



SO that 



m 
M 



where we have redefined W{x) again to include the scaling factor x~2. 

To begin, we open the divisor function d{n) and let n = nin2. To locate the variables rii 
and n2, we apply a dyadic partition of unity to both rii and n2 so that rii x Ni and n2 N2, 
where N1N2 x A^. Without loss of generality assume N2 > Ni. Let Bm,Ni,N2 be the sum 
B'jyj but with this partition of unity applied, and similarly for Am,Ni,N2{''t^)- 

We prove the following in Section IT2] 

Lemma 4.1. Suppose N2 > Ni, N1N2 N :$> M, MN < g^+e^ j^j^^^ 
(4.2) Bm,n„n, < <j (giiV, ^ + ) ^ N^q' 

The first bound is nontrivial for NM~^ > q^^^ while the second bound is nontrivial for 
g2+^ < Ni < NM~^q~'^ (note this upper bound is achieved automatically if ^ M^g"^ since 
Ni < y/N). The third bound is nontrivial for Ni < N^M-^q-". Note that if NM''^ < q, 
then there is a gap that has not been covered. In Section 14.31 we treat the case where Ni is 
relatively small, proving 



Lemma 4.2. Suppose N2> Ni, N1N2 N :$> M, MN < Then 

(4.3) Bm,n,,N2 < fiV^min (^^,^=] + m'^ N'"^ + mI N'! N, | q'. 

N sMs V MN I 



We briefiy describe how to deduce Theorem 13.21 from Lemmas 14.11 and 14.21 Basically we 
need to find the maximum value of the minimum of the above bounds as A^i varies between 1 
and VN. The case A^i = 1 gives < N'^M^q"" and A^i = a/ZV gives < (M^A^'i + g3A^"i)g^ 

1 1 3 3 

(which already dominates the Ni = 1 contribution). Note that M2N~2 + M^N'^Ni is 
dominated by NiM~2 provided M < A^a, which we may assume since otherwise (I3.15P is 
trivial. Hence the problem reduces to finding 



-- ■ i . / 92 



max (min(ciA^i ^,C2A^i)), where Ci = g* C2 = min , . . . , , 

i<m<VN^ \NlMiWMN 



The maximum is clearly 
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4.2. Proof of Lemma 14.11 We begin by noting that we may assume {n, d) = 1 since the 
contribution to Bm,n from q\n is 0(g~^^^) which is smaller than all the bounds in fl4.2l) as 
well as (14.31) . Define A' to be the same as A but with (n, d) = 1. 

Next we break up the sum over n2 into arithmetic progressions (mod d), getting 

(4.4) Av„.,.,(™) = -i-.E«/'') E E "'■(T)^(^)"'(i) 

d\q (ni,d)=l n2=rrmi (mod d) 

Using the same argument we used to remove V reduces the problem to estimating 
A'j[^^^ ^^(m), where A" has the same form as A' but with W{nin2/N) removed. Apply- 
ing Poisson summation to the summation over n2 gives 

-«,..^.M-^E^(« E f Ee(^)>^(^)5^(^). 

^ d\q {ni,d)=l h ^ ^ \i/ \ I ^/ 

where W is the Fourier transform of W . The term h = contributes to A'^^ N2 precisely 

N2 



^-(o)E--(^)E^.(,/.)«^/ 



ni ^ d\q 



and hence its contribution to B'^j ^ is 



< 7;—q' < q 

g2 



From the rapid decay of W, the terms with d = 1, h ^ contribute 0{q ) to Bm,Ni,N2- 
It remains to consider the terms with d = q, h ^ 0, namely 

where H = qN2^ . Due to rapid decay of the Fourier transform, we may assume \h\ < Hq^ 
(which implies N2 < g^~^^). Furthermore we may assume {hm,q) = 1 at no cost. 

It is possible to bound the sum over rii using the Weil bound for Kloosterman sums, 
however one can obtain better bounds on average over m and h. We glue together m and h 
to create a longer variable / = mh. We have 



E-^(^)-(iT 



(4.6) fl«J^V- 

where L x MH. Bounding this double sum is of independent interest. Friedlander and 
Iwaniec ( |FrI] . p. 337) state a bound for this sum which is nontrivial for A'^i and L rather 
short (e.g., L = q^^^ and Ni = g^/®"*"^), however we do not use their bound. We have 

Proposition 4.3. Let W be a smooth function with support in [1,2], and let 

(4.7) S{K,L;q):= J^' 



0<1<L 

IfK,L<t: then 

(4.8) ^(i^,L;g)« 



E-Mf)-(| 



Lq 



+£ 



L-2q4, + K'iL)q^ 
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Delaying the proof of Proposition 14.31 for a moment, we complete the proof of Lemma WTLl 
by inserting fl4.8p as well as the trivial bound into fl4.6p . 

In the course of proving Proposition 14.31 we need the following. 



Lemma 4.4. For any odd prime q, and x, y, z, G ¥q, let 
(4.9) 

Then we have the following bounds, 



T{x,y,z;q):= ^ ^ ^ e 

a (mod q)b (mod q)c (mod q) 



c{a — b)\ f ax + by -\- cz 



(4.10) 



T{x,y,z;q) < < 



q, z = 0,x^-y, 
g^ z = 0,x = -y, 



g, z^O,xy = 0, 

3 



x,y,z ^ 0. 

Actually the proof gives an exact formula for T in terms of Kloosterman sums. 

Proof of Lemma \4-4\ First note that the sum over c is — 1 unless a — b + z = 0, in which case 
it equals q — 1- Hence 



T{x,y,z;q) = (q-l) ^ ^ ^ |^axj^^ 



a~bz=~z (mod q) 



a—b^—z (mod q) 



ax + by 

q 



Solving for b and using a + z = a(l + za) gives 

ax + a(l + za)y 



T{x,y,z;q) = q e( 



E* v^* f ax + by 
1^ ^ 



a (mod q) 
a^—z 

where in case we interpret the condition a ^ —1 to be vacuous. If q\z we get 

T(x,?/,0;g) = qCq{x + y) - Cq{x)cq{y), 

where Cg(n) is the Ramanujan sum. To continue the calculation we impose the condition 
{z, q) = 1. Using a ^za — 'z and calculating the sum over a as a Kloosterman sum, we get 



T(x, y, z\ q) = qS{xz, -yz; q)e 



-xz + yz 

q 



q - Cg{x)cq{y). 



Finally, we obtain f l4.10p from case analysis, the Weil bound, and the easy computation 
Cg{n) = g — 1 if q\n, and — 1 if g f n. □ 

Proof of Proposition \4 ■ 3[ Applying Weil's bound for Kloosterman sums to the sum over k 
(by completing the sum) gives 

(4.11) S{K, L; q) < Lq^+' + KLq-^ < Lq^+' . 

Further gains can be obtained by using the sum over /. To begin, we smooth the sum S to 



S„(A-,L) := 



w 



Wo 
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where Wq is a nonnegative, Schwartz-class function satisfying ti?o(x) > 1 for < a; < 1. By 
positivity, S < Sq. Applying Cauchy's inequahty gives 

,4.12) s.iK,Lr « le* r r ^ (^) - (I) - (I) (l) ■ 

Now apply Poisson summation in each of /, ki, and k2 (mod q) to give that the right hand 
side above is 



From ( I4.10p and analysis, we get 



q3 LK^ ' ^ K 

4.3. An estimate for Ni small. In this section we prove Lemma 14.21 We continue with 
the estimation of f l4.5p . recalling that the error terms prior to arriving at (14.51) are acceptable 

i-T— 1 13 1 

for Lemma [4.2[ It is instructive to consider the crucial case M = q^ ^ N = q2 ^ Ni = q^ , 
N2 = q (whence H = 1). We begin with the simple reciprocity law 

f hmn{\ f—hmq\ / hm\ 



\ q J \ rii J \qni J 
where qq = l (mod Ui). Using e(^) = 1 + O(^) and = q^^ + 0{q^'^) gives 

Let Ri be the new sum above. Let g = {h,ni), whence 

c/(fc.»,,) = E<'('")»(^)>^(i)- 

Using the Mellin transform of W to write the sum over m in terms of the Estermann 
function gives 

(4.13) U(h,m,ni) = —[ M'D(s,0,^)W(s)ds. 

2m J(2) rii 

Next move s to —1, crossing a double pole at s = 1 which gives a 'main term' of size 
^ Mn^^q'^, which contributes ^ M2N~2q^ to -Ri, an acceptable error term. Apply- 
ing the functional equation of the Estermann function (Lemma 12. 7p gives U{h,m,ni) = 
U+{h, m, ni) + U-{h, m, rii) + 0{Mni^q^), where 

(4.14) U±{h,m,ni) = — [ M'n\-'^'D{1 - s,0,^)W±{s)ds, 
where 

W±{s) = W{s)2{27rf'-^T{1 - sfC±{ns), 



THE FOURTH MOMENT OF DIRICHLET L-FUNCTIONS 19 

and C_|. = 1, C_ = —cos. Expanding D into absolutely convergent Dirichlet series and 
letting s — > — s + 1 gives that 

U±{h, m,ni) = — > d{k)e y± , 

\ ni J \ni J 

where 

y±{x) = ±— / x-'W{l - s)2{27i)-'T{sfC±{7fs)ds. 

Note x^y^\x) <^a x"^ for any A > (by taking Re(s) large) and y±{x) = Cilogx + C2 + 

0{x^~^) for X small (by taking Re(s) = —1+e). Also y±\x) = cj + 0{x^~'^) for j = 1, 2, . . . 
and X small. 

Thus Ri = R+ + R_ + 0{M^N-^q'), where 
(4.15) 

where here Wq{x) = x~^W{x). The trivial bound gives 

(4.16) i?± < -F^q"- 

It now suffices to show 



which we proceed to prove. Let V{h, k) be given by 



v('a,^E^e(f).(-).„(^). 

Our next goal is to use the theory of exponential sums to prove that there is cancellation in 
the sum over rii, and bound the sums over h and k trivially. The presence of /i is a bit of a 
nuisance, so we again use the elementary reciprocity to law to write 

kqh\ f —kqnl\ ( kq 



rii J \ h J \hni 
By splitting rii into residue classes (mod h), we get 

We may assume gh <ti Hq'^ = qN2^q^, which is ^ Niq~^ since our results are trivial if 
N <^ g^^^ (the point is that there is room for summation over r). 
By partial summation it suffices to bound 

kq 



h(a + hr) 
id 

reproduce here (with slight changes in notation) 



gh — gh 

for C <^ We apply the convenient bound given by Corollary 8.5 of [IK], which we 
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Proposition 4.5. Suppose f{x) satisfies 

for k = 2,3 in the segment [Q, 2Q]. Then for 1 < Q' < Q < F we have 

J2 e(/(m)) < AF'eQh logSQ, 

Q<m<Q' 

where the implied constant is absolute. 
For our application, 

^^^^ h{a + hx)' N,h' ^ gh' 

and A is absolute. Then we get 



Using this bound gives 



^'■''^ "-^ (i) [m) (7) ' ■ 

which simplifies to give fl4.17p . as desired. This completes the proof of Lemma [4. 2[ 

5. A BINARY DIVISOR SUM 

When M and N are roughly the same size we treat Bm,n as a binary divisor sum. We 
loosely follow the presentation of Motohashi, Chapter 4 |M2] . Before beginning the treatment 
of the divisor sum, we first present some material to which we shall shortly refer. 

5.1. Some arithmetical sums. We require the computation of some arithmetical sums. 
Lemma 5.1. Let Re{s — A) > 1 and Re{s) > 1. Then if d is either 1 or a prime then 

(5.1) E ^-CWC(.-AM-E^->W- 

n=0 (mod d) bc=d 

Actually Lemma 15.11 holds for general d but we only need it for c? = 1 or prime. In this 
special case the proof is greatly simplified so we omit it. 

Lemma 5.2. If Re{a) < we have 

ci{n) 



(5.2) ^^^n) = ai-a)Y^'f^ 



-a 



Proof. This follows immediately from using the sum of divisors formula for the Ramanujan 
sum, that is Q(n) = Ylid\{in) dfi{l/d), and reversing orders of summation. □ 

Lemma 5.3. Let d be either 1 or prime, Re{s) > 1, and Re{X) > —1. Then 



/2+A dX(2 + X) V 

/ h (mod n=0 (mod d) \ J \ 



1 + 
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Proof. Evaluating the sum over h gives that the left hand side above equals 

1 ci(n) 



(5.4) E E 



p+x ■ 

n=0 (mod d) I 

Executing the sum over / using fl5.2p shows it equals 

(5 5) ^ V ^-^-"^"^^ 

C(2 + A) ^ • 

^ ^ ^ ri=0 (mod d) 

Finally using Lemma 15.11 and the fact that d is either 1 or a prime completes the proof. □ 

5.2. An approximate functional equation for the divisor function. 
Lemma 5.4. For any positive integer n and A G C, 

(5.6) -.(") = E 7^/^(7^) + « ^ 7i^^-<7i>' 

(5.7) fx{x)= [ x-'^Cil - X + w)^^dw, 

J{a) w 

a > |i?e(A)|, ci{n) is the Ramanujan sum, and G is as in Definition \2.5l 
Remarks. 

• The reason to expand the divisor function into such a series is that the exponential 
sum formula for the Ramanujan sum ci{n+f) will allow for the separation of variables 
n and / in a\{n + /). This is a simple alternative to the delta method [DPI] . 

• Motohashi used the formula (15. 2p to accomplish the separation of variables. However, 
(15. 2p does not hold for a in a neighborhood of the origin. 

• An approximate functional equation is essentially equivalent to a functional equation. 
The sum of divisors function ax{n) satisfies the functional equation 

axin) = n^a^xin). 

The function n^^ax{n) is perhaps more natural to study because it is invariant under 
A ^ —A. Of course, these appear as Fourier coefficients of Eisenstein series (12. lip . 

Proof. By a contour shift we have 

f ax-Un)n^^'^dw = ax{n)+ [ ax-Mn'^'^^dw. 

J{a) W W 

An application of the functional equation (Jx-w{n) = n^~^o'^x+w{n) and the change of vari- 
ables w —w gives 

axin)= [ ax-Mn'"/'^dw + n^ f a^x-^{n)rf''^^^dw. 

J(a) W W 

Inserting (15. 2 p into the above integral representation and reversing the order of summation 
and integration completes the proof. □ 
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5.3. Separation of variables. We begin the treatment of Bm,n by solving the congruence 
m = n (mod d) by letting n = m + f , where d\f, and / > 0. We have 

aa-(3im)a^^sim + f) f rri{m + f)\ f m m + f 



ma 



Our immediate goal is to separate the variables n and / both arithmetically and analyti- 
cally. We use the expansion of a^-s{n + /) into Ramanujan series given by Lemma 15.41 to 
arithmetically separate the variables n and /. Let Cm,n be the contribution from the first 
term of Lemma [5.41 and Cm,n be the second term. Then 

(5.8) Bm,n = Cm,n + Cm,n, CA/,Ar(a,/3,7, (5) = CM,Af(a,/3,(5, 7), 

where the formula for Cm n is 



dlci I 



d\q I f=0 (mod d) 



aa-f^{m)ci{m + f) ( m{m + f)\ [ m m + f\ f I 



EUa-pyiiijCiyni -r J ) / inyin -T J ) \ in in T J \ , 

— T— — — T— V7, — /' 



The function should not be confused with the variable of summation /; at any rate, this 
alphabetical accident shall clear up shortly. Before stating our preferred formula for Cm,n 
we set some notation. Define the Mellin pair 

We generally use the notation J^^ ^ to denote the vertical line of integration with Re(s) = Cg, 
and similarly for other subscripts. Let 

(5.10) H,{s, ui, U2, w) = ^^^}^^g{s)W{ui, U2)C{1 --f + 6 + w), 

sw 

(5.11) H{S, Ui, U2, V, W) = — ^— ■ ■ — Hi{s, Ui, U2, w). 

r(2 + 7 + S + M2- 2") 

We claim both Hi and H have rapid decay as any of the variables get large in the imag- 
inary direction. This is easy to see for Hi. For H we note that a crude use of Stirling's 
approximation shows that (for b and v in some fixed vertical strip) 

(5.12) '^^iV'^^ ^^^ + Hr^a + Hr, 

where A is any sufficiently large (depending on the fixed vertical strip) positive number, and 
for b and v avoiding the poles of the gamma function. For example, if \v\ < 2\b\ then the 
ratio of gammas has exponential decay in v, uniformly in b. In the opposite case just use the 
fact that it is bounded by a polynomial in v and b and multiply and divide by (1 + 1^1)^"^. 
The rapid decay of Hi in all variables overcomes this potential growth in b to show H has 
rapid decay in all variables. 
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Lemma 5.5. With Cg = 2, = = e, and c^^ = = 0, we have 



^ d\q I h (mod /=0 (mod d) J ^ ^ (Cs) J (c^) 

/ / / — - — ^ — Dl'^ + a + s + Ui + v,a — P, —)H{s,Ui,U2,v,w)dvduidu2dwds. 

J{c^,)Jic^,)J{c,) " ^l"" I 

Proof. Inserting the exponential sum formula for C;(m + /) into fl5.9p gives 

^ d\q I /i (mod /=0 (mod d) ^ 



Next we separate the variables m and / analytically by taking Mellin transforms. We have 

(5.14) ^^y{^^]U^)w(^r^^ 



I ^ : / / / / 1 —du2duidwds. 

\27iiJ J^,^^) l^m^+^^{m + /)2+7+-+«2-2 

To separate m and / in (m + /)~^~'T~^~"2+f ^j-^g following formula (17.43.21 of |GR] ) : 

(5.15) (1 + -)-' = ^/ " 

valid for < c„ < Re(6). Recall fl5.12p for the convergence. 
Thus fl5.14p equals 

(5.16) I I [ [ [ [ [ - — 1 — - — - — — ^-dvduidu2dwds. 

Inserting (I5.16P into Cm,n, we obtain 

^ d\q I h (mod /=0 (mod d) J m 



m2" 



I \^ [ f f f f q^'M'''N''^H(s,Ui,U2,v,w) , , , , , 

-dvauidu2dwas. 



2vr^; Jic.) J(e.) i(c„,) i(c„,) i(c„) /^+"2-.--/»^s+..+. 

The sum over m converges absolutely on the contours of integration as stated in Lemma [5751 
Writing this sum over m in terms of the Estermann D-function finishes the proof. □ 

Our goal is to develop the asymptotics of Cm,n by moving the lines of integration past the 
poles of the D-function (which will contribute the main terms) and to analyze the remainder 
term using estimates for sums of Kloosterman sums using the Kuznetsov formula. 

To begin, we move the s-line of integration to Cg = e, passing the two poles of D at 
| + s + Mi + f = l — a and ^ + s + Ui + v = 1-/3. Let Em,n{c(, (3, 7, 6) be the integral on the 
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new line of integration, and let Pm,n = Pm,n{c(, (3,1,^) be the contribution of the former 
pole; by symmetry considerations, the latter pole is PAf^Ar(/3, a, 7, 5). Thus we have 

(5.18) CM,7v(a, (3, 7, S) = Pm,n{o:, (3, 7, S) + Pm,n{P, a, 7, ^) + Em,n{o:, 13, 7, 5). 

We bound Em,n in Section [HI proving Theorem 13.31 and continue with our calculation of 

Pm,n in Section O 

6. The main terms 

The partition of unity is an obstruction in the computation of these main terms. It turns 
out to be easier to sum over M and before doing finer analysis of the main terms. 

Lemma 6.1. For N ^ M, we have 

(6.1) Pm,n <^ N^'^M'^q^ 

We prove Lemma 16.11 in Section 16.21 

6.1. Recomposition. Define P := Pm.at- 

Lemma 6.2. We have 

^G{v)G{w)g{v) 



vw 



(6.2) P(a,/3,7,5) = C(l-a + /3)(^^^ J^^ ^ q'^ 
V{\-a-v)T{a + ^ + 2v-^) + ^ + - f )C(1 + /3 + (5 + 2t; + f ) 

r(i + 7 + ^-f) ^(1-^+^+-) c(2-a + /3-7 + ^ + ^) 

+ 0(g-i+^). 

The four main terms from B'q are obtained by taking the main terms of Bq and switching 
a and 7, and (3 and 5. One term is 



r(i - 7 - v)T{a + 7 + 2i; - f ) C(a + 7 + 2i; - f )C(1 + f3 + 6 + 2v + ^) ^ , 
= = = = — nun HI 

rC^ + a + v-f) C(2-a + /3-7 + 5 + ^) 

and similar formulas hold for the other terms. 

Proof. We first show how to evaluate the sum over M and N. As a warmup problem, we 
show that for a 'nice' function F, 

( 7T— ) / / F{ui,U2)WM,N{ui,U2)dUidu2 = F{0,0), 

^-^ \ 2m J J(,. ) J -s 

For a proof, let / be the inverse Mellin transform of F and begin with the Mellin convolution 
formula to write the left hand side above as 

t^Jo Jo X y Jo Jo X y 

The point is that on the level of the Mellin transform, we divide out by W{0, 0) and evaluate 
everything at Ui = U2 = 0. 

The computation of P is of this form but with the extra step of computing the residue of 
the Estermann function at| + s + f = l — a. We make a brief diversion to justify this step. 
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Note that Pm,n is the difference of two integrals of the form (15.131) with the only difference 
between the two being the placement of the line of integration c^. Applying the changes of 
variable s — )■ s — mi and w ^ w — 2ui + 2u2 transforms f l5.13p into 

^ d\q I ?t (mod /=0 (mod d) ■> ^ ^ -J {Cs) J (c^) J (c^^) 

/ / j5 D(^+a + s+v,a — B,—)H2(s,ui,U2,v,w)dvduidu2dwds,dean 

where 

, r{v)T{l + j + s-v-^)G{s-Ui)G{w-2ui + 2u2) 

H2{s,Ui,U2,V,W) = ■ 



r(| + 7 + s-f) {s-ui){w -2ui + 2u2) 

g{s - ui)W{ui, M2)C(1 - ^ + 6 + w-2ui + 2u2). 
Hence after summing over M and in ( 16. 4p we get 

^ d\q I h (mod /=0 (mod d) J 

{^~] III .s-!-H„„ -P(| + a + g + j)H3{s,v,w)dvdwds, 

where is the following function with rapid decay in all variables (recall (I5.12p ): 

X ^Wih + l + s-v-^)G{s)G{w) , 

6.6 H^{s, V, w) = ' \ ^ ^ ' ^ ' g{s)C{l -l + 5 + w). 

r(i + 7 + s--|) sw 

Thus P(a, /3, 7, 5) is the residue of the integrand in (16. 5 p at s + f = | — a. 

Using Lemma 12.71 we find that the residue of the Estermann function at this point is 
/-i+"-''C(l - a + and hence 

^^^W^5:«/'')Ei5^E* E 

^ d\q I ft (mod /=0 (mod d) ^ ^ 

I I / / — ; \ — n — s — H^i^ — a — ihv,w)dvdw. 

Apply the change of variables v^^ — a — vto get 

^-^^^^E«/'')Ef^E* e if) 

^ d\q I /i (mod f =0 (mod d) ^ 



d\q I /i (mod /=0 (mod d) 

Hsiv, ^ — a — V, w)dvdw. 



.2vr^; 7(c„)i(c„)r+^+'''-t/ 

Note that the ratio of gamma functions in above is 

r{l-a-v)r{a + -f + 2v-f) 

r(| + 7 + ^-f) ■ 
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Since — a) = (recall Definition I2.5p . we may move to 1 without encountering any 
poles. The summation over / converges absolutely, and we may execute the summations 
over /, h, and / using Lemma [5.31 to obtain that 

(6.8) P = ~ ^ + <i^id)Md) (^) ' I I q'''H,{v,l-a-v,w) 



d\q 

C(a + 7 + 2t; - f )C(1 + f3 + 5 + 2v + ^) 



rf°+^+2''-2((2-a + /3-7 + 5 + w;) ^ ' 
Now we move the f-line of integration back to = e. The pole ata + 7 + 2i; — ^ = 1 gives 



^ d\q 



(1 + ^-1+- /3+7-<5— _ ^-2+a-/3+7-5-«,) ^^(1 _ « - 7 + i(-a + 7 - -^),w)dw, 

which is seen to be 0(g~^^^) due to cancellation in the arithmetical sum over d. Using trivial 
estimations shows that 

P = C(l-a + /3)f-^V / / 



.2vrv j(c„)j(c„) 

—H-i[v, j: — a — v,w) + 0[q ^ ' 



C(2 - a + /3 - 7 + (5 + w) 
Observe that 

vw r(^ + 7 + f-|) 

which completes the proof. □ 

6.2. An upper bound for Pm,n- In this section we prove Lemma [6.11 A minor variation 
of some calculations in the proof of Lemma 16.21 reduces the problem of bounding Pm,n to 
bounding 

C(l - a + /3) f-^) [ [ [ [ g-"-^-"^-"'+^M"^A"W(Mi,M2) 
\2mJ J (c^) J (c,) J (cu^) J (c^^) 

C{a + -f + Ui+U2 + 2v- f )C(1 + (3 + S + ui + U2 + 2v + f ) G{v)G{w) , . 



C(2 - a + /3 - 7 + 5 + to) vw 

C,{\ — 7 + (5 + vS)du2du\dvdw 



r(i - a - ui - v)Y{a + 7 + Ui + Ms + 2t; - f 



2 I / I - I "-z 2 I 



r(i + 7 + t; + U2- 

where = c„2 =0; c^, = = e. Then move c„j to ^ — 26 followed by c„2 — ?■ — |, and 
bound the integrand trivially with absolute values to finish the proof. 

7. The dual off-diagonal terms 

To properly manipulate the main terms obtained in Section |6l it is necessary to combine 
them with the corresponding main terms of Aq. The computation is similar to that of Bq 
but there are some differences. 
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As before, we require different methods of estimation depending on how close m and n are. 
We apply the same partition of unity as before and write Aq = J2m n ^m,n accordingly. 
Clearly Am,n satisfies the bound f l3.15p . 

7.1. Separation of variables. We follow the methods of Section |5l Write f = m + n so 
Am,n = — ^ Yl 't>{d)^^{(l/ d) Yl 

^ d\q /=0(mod d) 



p{m)a^_s{f - m) f m{f - m)\ f m f-m 



We use the expansion of <J-y-s{f — m) into Ramanujan series (Lemma 15.41) to separate the 
variables. Write Am,n = Cjnr + Cjjj^ to correspond to the first and second parts of the 
approximate functional equation. Then 



d\q I /i (mod i) /=0 (mod d) ^ 



^ ua-f}{m)e{^) r m{f - m) \ f m f-m 



'MM 

^ d\q I h{modl)f=i 

/>0 

0<m<f \J / ^ ' ^ 

Using the same methods as in the computation of Cm,n, we obtain 

^ d\q I /i (mod /) /=0 (mod d) m 

f 1 Y r r r r q^''M'''N''^Hi{s,ui,u2,w 



,11111 1 „ dsduidu2dw, 

V27rV J(c^,) J^c^^) J(c^^) J^c^) /-m«+«i(/ - m)2+T+-+«2-2 

where Hi is given by f lS.lOp . 

To separate the variables m and / in (/ — m)^^^'^^^^'^^^^'^^ we use the following formula: 

1 f T(v)T(l-b) , 0<x<l 
-X dv — 



2m J^^^)T{l-b + v) [0, l<x. 

valid for < c^, Re(6) < 1, from 17.43.22 of [GRj . The integration converges absolutely 
provided Re(6) < using Stirling's approximation, namely 

^'f'\-''' «(l + |I.nMir'>)-', 

r(i — b + V) 

where the implied constant depends polynomially on b (the dependence on b could be ex- 
plicitly stated but it turns out not to be relevant in this work). Thus we have 

^ d\q I h (mod I) f=0 (mod d) J m 

1 Y f f f f r g2sjvf«iAr«2ms,Mi,M2,w,w) , , , , , 

-dvduidu2dwds, 



2vr^; 4.) i(c„) J(c.,) Jic.,) Jic.) fs+u,-.-^i^^s+u,+v 
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where 

H{s,Ui,U2,V,W) = —J ■ ——Hi{s,Ui,U2,W). 

- -f-S-U2 + V + f) 

Note that this expression is identical to fl5.17p except that the ratio of gamma factors arising 
from the separation of variables is different. We essentially follow the same outline used to 
compute Cm^n but some arguments must be altered due to changes in location of the poles 
of these gamma functions. We initially take Cs = |, = = 0, = e, and = 2 + e so 
that all the integrals and the sum over m converge absolutely. Then we write the sum over 
m in terms of the Estermann D-function and move Cs to e, thereby crossing the two poles of 
D. Let -Ejvnv be the integral along the new line of integration, and Pju^{a, (3,'y,S) be the 
contribution of the pole at ^ + s + ui + v = 1. We compute -E^nv Section [9] and proceed 
with the computation of the main terms. 

7.2. The main terms. The same argument used in the proof of Lemma [6.21 shows how we 
may execute the summation over M and A^. Furthermore, with the same method of proof 
as Lemma 16.11 we have 

(7.1) P^«M^iV-ig^ 

Let P(a,/3,7,(5) := J2m n ^Mj^i^^ f^^'^^^)- "^^^ analog of Lemma 16.21 is 
Lemma 7.1. We have 

(7.2) Pia,P.-,.S) = ((l-a + P){±-)' f f ,-°-^^f ^'"'^.'"''^'"' cd - 7 + ^ + ^) 



(e)Jie) 



vw 



r(i - g - t;)r(i -j-v + f)C{a + j + 2v- f )C(1 + (3 + 6 + 2v + 
r(l-a-7-2^; + f) (^(^2 - a + /3 - ^ + 6 + w) 

+ 0(g-i+^). 

Note that the pole of ({a + •y + 2v — ^) is cancelled by one of the gamma factors. 
Proof. To begin, P is given by the residue of the following integral at s + f = | — a: 

^ d\q I /i (mod /=0 (mod d) J 

(^■] Iff ..-!-H„„ -P(| + a + s + v,a- j)H^{s, v, w)dvdwds, 

where 

T7 . ^ r(t;)r(i - 7 - s + I) G{s)G{w) , , 

^3 V, w) = -^^^ ' ^ ' ^ ' 9{s)a^ -y + 5 + w). 

Now we compute /3, 7, 5). Computing the residue and changing variables v \—a—v 
gives the following analog of (16. 



^ d\q I /i (mod /=0 (mod d) ^ ^ 

I I / / — —77^ — w — Hr,(v, ^ — a — v,w)dvdw, 
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where note the ratio of gamma factors in i^s is 



r(l - a - 7 - 2^; + f ) 

Now we take = and move to \ + e. The pole of r(| — a — v) is cancelled by a zero 
of G, and the sum over / converges absolutely, so we may borrow the same computations 
used in the proof of Lemma 16.21 to get fl7.2l) . □ 

8. Assembling the main terms 

We now begin to assemble the various main terms to form a nicer expression. We briefly 
recall how the main terms decompose. Recall (13. ip . (13. 2p . and (13. 4p . In general, (13. 2p shows 
how to derive results for the "second part" of the approximate functional equation. Recall 
the combinatorial dissection (13.110 . as well as the computation oi Ao from Lemma 13.11 
Then we wrote Bq = J2m n ^m,n and recall (15. Sp . (IS.lSp . and Lemma I6l2l This gives Bq 
as the sum of four main terms, and we similarly express Bq^o, /3, 7, 6) = Boipf, ^, a, P). The 
computations are similar for the dual terms as in Section [71 



8.1. Combining Aq and Aq. To begin, let Q = Q{a, ^,^,6) be the sum of (16. 2p . (16. 3p . 
and (17. 2p . From the above discussion, the total contribution of main terms from Bq, Bq, 



and Aq is then 



•1) 



Q(a, 13, 7, S) + Q{f3, a, 7, S) + Q{a, (3, S, 7) + Q{f3, a, S, 7). 



This grouping of the main terms is suggested by unpublished work of Hughes on the Riemann 
zeta function [Hi. 



Lemma 8.1. We have 



1.2) Q(a,/3,7,^) 



C(l-a + /3)C(l-7 + ^) 
C(2-a + /3-7 + 5) 

C(l-a-7-2s)C(l + /3 + 5 + 2s) 



r 



7— s 



r+Q + S 



I-7 + S 



■ds 



Proof. By its definition, 
(8.3) Q 



1 

27Ti 



-o— 7+- 



Giv)Giw)giv] 



vw 



7 



R{v, w) 
+ 2v- 



':)({! + P + 6 + 2v 



2 ^ 



where 

R{v, w) 



+ Cil-a + (3 + w)C{l-i + S) 



C(2 - a + /3 - 7 + (5 + w;) 
r(|-a-^)r(a + 7 + 2i;-f) 
r(i-7-t;)r(a + 7 + 2t;-f 



dwdv. 



r(i + a + t;-f) 



T{1- a--f -2v + 
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We moved the line of integration to \ without encountering any poles. Now move the line 
of integration over w to —1 + e, crossing poles at w = 0, and ^ = — | + 7 + f only (recall 
Definition I2.5p . The contribution from the pole at w = — 1 + 27 + 2f is 

It is clear that this integral is bounded by 0{q^^/^^^)^ as can be seen by moving f to e and 
capturing the pole at 3w ~ |. 
The pole at w = gives 

2injfiy V C(2-a + /3-7 + i) 

Now apply the functional equation to C(a + 7 + 2v) to express this as 

1 f f q \ -"-7 G{v)g{v)R{v, 0) r( ^"""^""'^ ) ^(1 _ ^ _ ^ _ 2^)C(1 + (3 + 6 + 2v)_^^ 



27rz7(i)V7ry n-2-'-v r(^^±^) ({2 - a + 13 - ^ + S) 

Note 

Riv, 0) = C(l - a + /3)C(1 - 7 + 5)r«,^(^^), 

where 



^.5) Ta,^{v) 



r(i - a - w)r(a + 7 + 2t;) r(i - 7 - v)r{a + -f + 2v) 
T{^ + -f + v) ^ T{l + a + v) 



T{l-a--f-2v) 
To finish the proof of Lemma 18.11 we use the crucial identity 

' ' - r (i-^) r (i±|±i) r (i±i±i) ' 

which is deduced from the following Lemma with a = ^ — a — v and 6=| — 7 — v. □ 
Lemma 8.2. For any a, 6 G C, a,b,a + b ^ Z, 

(8 7^ r(a)r(i-a-&) ^ r(&)r(i - a - 5) ^ r(a)r(6) _^xr(i^) r(f) r(|) 



r(i-6) r(i-a) T{a + b) T {^) r{^)r{^)' 

Proof of Lemma \8.^ Using 



the right hand side of (18.71) equals 



l-s 



(8.8) 2- 



cos(f )cos(f )r(a)r(6) 



cos 



:^^)r(a + 6) 
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Using a series of standard gamma function and trigonometric identities, the left hand side is 
T{a)T{b) fT{a + b)T{l- a-b) , T{a + b)T{l - a - b) ^ ^ 



T{a + b) \ T{b)T{l-b) T{a)T{l-a) 
r(a)r(6) / sin(7ra) + sin(7r6) + sin(7r(a + b)) 



2 



r(a + 6) \ sin(7r(a + 6)) 

r(a)r(6) sin(f (a + b)) cos(f (a - b)) + sin(f (a + b)) cos(f (a + 6)) 
r(a + 6) sin(7r(a + 6)) 

= 4 



T{a)r{b) sin(f (g + b)) cos(f ) cos(f ) 
r(a + b) sin(7r(a + b)) ' 

which simphfies to give (18. 8p from the sine double angle formula. □ 

8.2. Combining Q and Q^. In this section we calculate the terms corresponding to Q that 
create the main terms of M_i. Recall that these are obtained from Q by switching the signs 
of a, /3, 7, 6 and multiplying by Xa^^^'y^s- We combine Q{a, (3, 7, 6) and Q-{f3, a, 6, 7). 

Lemma 8.3. We have 

(8.9) Q{a, (3, 7, 6) + Q_ (/3, a, 6, 7) = U{a, (3, 7, S) + 0(g-i/3+^), 

where 

C(l - a + /3)C(1 - « - 7)C(1 + /3 + 5)C(1 - 7 + ^) 



(8.10) [/(a,/3,7,<5) = X, 



",7 



C(2 - a + /3 - 7 + 5) 



Note that f/(a,/3,7,5) is one of the middle four terms in Conjecture 12.11 Hence, adding 
(18. 9p according to (18. ip to the contribution of the diagonal terms given by Lemma 13.11 will 
form the quantity on the right hand side of Conjecture 12.11 

Proof. It is a matter of bookkeeping to modify ( 18. 2p to see 



C(l - a - 7 + 2s)C(l + /3 + 5 - 2s) ^ r^-l+s\ ^ } i-s+/\ 9^c.,-p,-,,-5{s)ds + 0{q-'/'+'^ 



(8.11) Q_(/3,a,5,7) = X„,^,,,5^^-^ ^ ±Tr^(-) 7^ 

C(2 — a + P — 7 + d) Vvr/ 27rz ^(i) s 



Let / be the main term in (18. 2p . and /_ be the main term in (18. lip . We work with / by 
moving the line of integration to — |, passing a pole at s = giving a residue which we easily 
compute to be f/(a, /3, 7, 5). Write I = U{a,P,'~f,6) + /', where I' is the new integral. We 
now show /' = — /_, which will complete the proof. 
Apply the change of variables s — )■ — s to give 

j,^ C(l-a + /3)C(l-7 + 5) 



C(2 - a + /3 - 7 + 5) 



7r 



G{s) 



— C(l - a - 7 + 2s)C(l + I3 + S- 2s) },./^ \/ t^-'^ g^,p,,,s{-s)ds. 
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Now we claim that 



q \ -"-7 V 2 7 V 2 



71" ^fl'a,/3,7,<5(-'5) 



which imphes /' = — /_ . Each of these terms is a product of terms each depending on exactly 
one of a, /3, 7, or 5, as well as a certain power of tt not depending on the shifts, so it suffices 
to check this identity at each such factor. The cases a and 7 are the same and so are /3 and 
5. The case of a follows from 



Y I izf^lf ^ Y ( Y ( ir"+'' 



Ij-a. — s\ _ /i4-fv\ ^2 ' _ /i 



2 



and the analogous formula for /3 is 



A-(| + «(i) 



/9 



We also check that the power of vr is the same on both sides. □ 

8.3. A note on odd characters. In this work we concentrate almost exclusively on the 
even characters in the proof of Theorem 11.11 The contribution of the odd characters carries 
through in the same way with slight changes. The only differences between the odd and even 
characters is that the X factors are different, and for the odd characters, the 'dual' terms are 
subtracted rather than added. The estimations of the error terms carry through as before; 
the only differences arise in the calculation of the main terms. Some thought shows that the 
evaluation of the main terms for the odd characters diverges from that of the even characters 
starting with the analog of (18. 5p : for the odd character case the third term is subtracted 
rather than added. The only essential difference is the use of the following gamma identity 
instead of Lemma |83 



Lemma 8.4. For any a,b ^ C, a, 6, a + 6 ^ Z, 

r{ani-a-b) mni-a-b) Tia)m _ , r(i^)r(i±^)r(i±^) 
^ ' r{i-b) ^ r(i-a) r{a + b) r(^) r(^)r(^) 

The proof is similar to that of Lemma 18.21 so we omit the details. 



9. Treating the error terms 

In this section we prove Theorem 13. 31 Since the forthcoming estimations become techncial, 
it may be helpful to know that the quality of the error term in Theorem 13.31 can be predicted 
by careful scrutiny of the formulas in Theorems 3 and 4 of [Ml] . 

In order to clean the upcoming formulas, we set all the shift parameters a,P, etc., equal 
to 0. All the arguments can easily be generalized to handle sufficiently small non-zero values 
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without degrading the results. Recall that Em,n is given by the right hand side of (I5.13P 
with contours of integration defined by 

(9.1) Cg Cyj £y C^(j c^j2 0. 

9.1. Reduction to Kloosterman sums. Apply the functional equation fl2.16p to D to 

obtain Em,n = + E_, where 

^ d\q I h (mod I) f=0 {mod d) J 

1 Y /• r r r r q^^M'^^N'^^ 

-2^; hcs) he.) he.,) hcu,) he.) Z^+'^-'^-f /2^+2«,+2.+. 

±h 
1 

and where 

(9.3) H±{s, Ml, U2, V, w) = 2(27r)-^+2^+2"i+2^r(| - s - m - vfH{s, Mi, Ms, v, w)S±, 
with 

(9.4) 5+ = l, S. =sm{'K{s + ui + v)). 

Here H± has rapid decay in all variables since H does and since the exponential decay of 
the gamma factors cancels the exponential growth of S^. 

Now move c^^ to —1 and expand D into absolutely convergent Dirichlet series and execute 
the sum over h in terms of Kloosterman sums to get 

5 



— s — Ml — t;, 0, -^)H±{s, Ml, M2, M, w)dvduidu2dwds, 



-H±{s, Ml, M2, M, w)dvduidu2dwds . 



As an aside, we mention that using only the Weil bound for Kloosterman sums, we obtain 
the bound 

Em,n <^q~^^'M-^N-i, 
which can be seen easily by taking 

Using MN ^ shows 

Em,n < M-^iVtg^ 

which is nontrivial only for N <^ M ^ g^"*""^. This is insufficient to combine with 

Theorem 13.21 to cover all ranges; for example, M = g^/^, N = q^l'^ would not be covered. 

To do better we shall obtain additional savings coming from cancellation in the sum of 
Kloosterman sums by the use of the Kuznetsov formula. 
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9.2. Preparation for application of Kuznetsov. Let 

'^^''^ - \2m) 4) 4.^) /2.+«,+..^f li^^j 

if± (s, Ui,U2, V, w)dvduidu2dwds , 

with contours of integration given by f l9.ip . Then 

^ d|g /=0{modd) m J '"' I ^ 

By taking = | — 2e, we get r±{x) <^ x^~'^. By taking c^^ = —A with A large, we see that 
r'^^^x) ^j^B (1 + x)^"^ for j = 0, 1, 2, . . . and any B > 0. These conditions are sufficient for 
the application of the Kuznetsov formula, Theorem 12.61 

The next step is to apply the Kuznetsov formula to the sum over I. We write E± = 
Em±+Ec±+Eh± to correspond to the Maass forms, the Eisenstein series, and the holomorphic 
forms (of course Eh^ = 0). We show how to estimate Em± and Ec±, since Eh+ is smaller 
and easier to handle (for example, see Section 5 of |M1] ). 

9.3. Integral transforms. At this point we manipulate the various integral transforms of 
r± that we require for the application of the Kuznetsov formula. We require Mr_^(t), Kr-{t), 
and Nr^{k), in the notation of jlKj, Theorems 16.5 and 16.6. We have 



/ 1 Y /• r r r r q^'M''^N''^H+{s,ui,u2,v,w) 

V2^; i(c.) i(c.) i(c„,) 4.,) 4.) (47r)2^+2«,+2.+»J2.+«,+..^f 

TT? dr 

{J2^t{x) - J-2it{x)) X^'+^'''+^''+^—dvdUidU2dwds. 



sinh 27rt ./n x 

Now use the formula 

^dx r(^ 



u--s±2 \ ' 



X r(^ 2 / 

vahd for | > Re(s) > -Re(z/) (see 6.561.14) of [GRj) to see that 



oo 



{J2it{x) - J^2it{x))x = 2 



X \^r(i-f + it)) r(i-f-2t) 

We simplify this expression using 
Lemma 9.1. 

T(a + ir) r(a — ir) 2i 

(9.6) — = smh(7rr) cos(7ra)r(a + ir)T[a — ir) 

V[l — a + ir) r(l — a — ir) vr 

We omit the proof since it is easy and standard. Thus 

7rz2^-l ( r(|+it) r(|-it)) \ 2^-1 „„„/^^N-n/^ , ..N-n/^ 



■cos( — )T{-+it)T{- - it). 



sinh27rt \^r(l - f + it)) T{l-^-it)) coshvrt ' 2 ' '2 ' ^2 
Hence 

Mr+it) = I I [ ... [ — w — +('^' ""2' dvduidu2dwds, 

' \2niJ J J ps+m+u2m^ cosh(7rt) 
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where 

ui, U2, V, w; t) = cos(7r(s + ui + v+ + ui + v+ — + tt)r{s + ui + v + — - it) 

and where c* is meant to account for bounded factors hke powers of 2, tt, etc. that do not 
have any effect on the convergence of the integrals. 
Similarly, we compute 

, , f I Y f f f f f q^'M^^N^'H_{s,u^,U2,v,w) 

Kr-{t) = \ — 

2s+2ui+2v+w ""^ 



2 / K2u{x)x'''+^'''+^''+'"—dvduidu2dwds. 



X 



This time we use the formula (see 6.561.16 of [GRJ) 



,_^_,,dx f s + 2it\ „ f s — 2it 

2it' 



Kou(x)x'— = r~'T 



X V 2 

valid for Re(s) > \Re{2it)\ = 0. Hence 

/ 1 V /" f q^'M'^^N'^'' - 

(9.7) Kr-{t) = -— ■■■ -z;-, — ; wH_{s,Ui,U2,v,w;t)dvduidu2dwds, 

\27n J J J 

where 

H_{s, u,v,w;t) = T{s+ui+v + —+it)T{s+ui+v + ——it)T(^ — s—ui—vysm{n{s+ui+v)) 

9.4. Maass forms. This section is devoted to proving 
Proposition 9.2. We have 

(9.8) « , 

We treat the opposite sign case E^- only since the case of Em+ is similar, and easier. 
Proof. We have 



which upon using (19. 7p is 



d{m) |p^.(l)|2A^.(m)A,(/) 



(9-9) E-- = J-X^idM'l) E E^E 

(p [q] , d ^ m'ii'i coshyrK 

^ d\q /=0 mod d m J ] 



... w-H_(s,Ui,U2,v,w; Kj) cosh.(7iK^)dvduidu2dwds. 

\2m J J J ps+ui+u2-ui2- 
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Let Ek be the same expression as (19.91) but with the spectral parameter kj restricted to the 

sans the su 

ax{n)\j{n) Lj{s)Lj{s — A) 



dyadic segment K < Kj < 2K. Taking = | and = | means the sums over m and / 
converge absolutely. Now we use the following variation on (13. 5p 



E 



«• ((2s - A) 

and 

E ^^.iE^)(E^ 

/=0 (mod g) •' (/,g)=l ">0 



to obtain 

(h*[q) ^-^ , g( coshvr/t,- \Zt:i J j(r\J(r \ 



r r C c?^ M^i Ar"2 



\ H_{s, Ml, ti2, W; Kj) cosh(7rKj) 



Ai((i) ) — — ^ — '-^ —dvduiduodwds. 

To estimate this term we initially move the lines of integration so that 

(9.10) Cg c-y C|^^ c, c ~t~ ^5 

where 3e < c < | — 3e, passing no poles in this process by close examination of the form of 

H_. Due to the rapid decay of we may truncate the integrals so that Im(s), Im(Mi), 
Im(M2), Im(f), Im(w) ^ {qKy with a negligible error (say, size <^ {qK)~^^^^) . The issue at 
hand is the dependence on K. Using Stirling's approximation, we see that 

(9.11) cosh(7rfi:j)r(s + Ui + v + ^ + iKj)T{s + Mi + w + f - inj) < q'K-^^^". 
Then we have 

(9.12) £,«,-lfeATA-«'(^)"'max| ^2 M^X^i,^L,(L + s,)LAi + s.) 

where the maximum over Si and S2 is over the rectangles < Re(sj) < e and |Im(sj)| ^ 
(gK)^^ = l,2. 

The spectral sum may be estimated using \j{q) ^ q^^^ and the fourth moment bound 

(9.13) V Mllt|i:,.(i + ,)|4«;^2+. 

^ ' ^ COshvTK,' -^^2 ^1 

K<Kj<2K ■> 

for Im(s) <^ K^, Re(s) > 0, which follows from the large sieve inequality for Maass forms 
|Iw2j . For a proof, see Theorem 3.4 in Motohashi's book |M2] (actually Motohashi had s = 
but this is a minor change). 

First suppose that K 3> M~^N^q^. Observation of the location of poles H_ shows that 
we can take c = A large without crossing any poles. For such K we thus obtain 

(9.14) Ek < (gir)-^°°°. 
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For the complementary range, taking c = ^ — 3e completes the proof with the bound 

(9.15) « g-^+^i^^ ' (gi^r « (^) ' ■ □ 

The estimate of Em+ is even easier than that for Em~ because cosh(7rt)Mr+()f:) has expo- 
nential decay as t — )■ oo. 

9.5. The continuous spectrum. In this section we prove 
Proposition 9.3. We have 

(9.16) i?e± « (^^ 

The exponent j depends on an estimate for the 6th moment of the Riemann zeta function. 
As in the previous section, we only show the full details for E^.- since is similar and even 
easier to handle. 

Proof. Our starting point is 
which after using (19. 7p becomes 

cr2it{rn)a2it{f) 



„.ifl7r7_^ (m/)^*|C(l + 2it)|2 

d\q /=0(modd) m J J oo \ J ) \^\ l\ 

-^H_{s,ui,U2,v,w;t)cosh.{'7it)dvduidu2dwdsdt. 



2vri J J J ps+ux+u2jy^2 

Now move to | and to 1 + e and execute the summations over m and / in terms of a 
product of zeta functions to get 

F rid r (1-^' 

d\q 



(Cs) J(Cw) J(Cu2) J {Cu^) J{Cv) 



^2s+u^+u2 ^-(^' ^5 i^) cosh(7rt) |^(-^^2it)|2 dvduidu2dwdsdt, 

where 

Z(s, Ml, U2, tt;, t) = C(| + f + ^t)'C(| + f - ^i)' 

C(| + 2s + Ui + U2 + it)Cil + 2s + Ui + U2 - it)Ad{s, Mi, U2, w), 

where Aii{s,ui,U2,w) is bounded by d^ and holomorphic for Re(s), Re('Ui), Re(M2), Re(w) > 
—6 (use Lemma [5TT] to get an exact expression). As in the previous section, take the contours 
according to (I9.10p where Se: < c < | — 3e is at our disposal. Although we crossed various 
poles of the zeta function, since they are all at height t (roughly speaking), the decay of the 
test functions shows that these terms have rapid decay in t-aspect, and it is not difficult to 
bound the contribution of these terms by ©(g"^*^*^*^). 
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Again, we may truncate all the integrals except the one over t at height {qty with negligible 
error. The crucial issue is convergence in t- aspect. Suppose we have a bound on the 6th 
moment of the Riemann zeta function of the form 



(9.17) y 



1 /"^^ 

The currently best-known result has 6' = 1/A using Holder's inequality with 

l>2T p2T 

Jt Jt 

This bound on the 12th moment was proved by Heath-Brown |H-B3] . 

Using (19. lip shows that c = —9' — e is sufficient for convergence in t-aspect, and gives the 
desired bound. If 6'' = (a standard conjecture) it would give the bound Ec± ^ g'a"*"^. □ 



9.6. The dual terms. In this section we sketch how to treat Ej-fj^. The basic outline is the 
same as for Em,n but some convergence issues are slightly more delicate and the arguments 
must be modified in some places. 

Using the same computations as in Section 19.11 we get that the analog of (19. 5p is 

1 y r r r r r q^sj^u^^u^^s+u^+v 

hcs) he.) he.,) hcu,) he.) /^+— f /2s+2«,+2.+. 

H±{s, Ui,U2, V, w)dvduidu2dwds , 

where H± is given by an expression identical to (19. 3p except H is replaced by H. More 
explicitly, the ratio of gamma factors 

(9.18) r(.)r(i + . + n.-.-f) 



r(i + s + u2-f) 

appearing in the definition of H± is replaced by 



(9.19) 



Furthermore, the term S{m, /; /) is matched with S- and S{m, — /; /) is matched with 5"+. 
The primary difference between (I9.18P and (19.191) is that the former expression has exponen- 
tial decay in f- aspect, and the latter has at best polynomial decay. 

Because the opposite sign case is matched with 5"+, the convergence in f- aspect is assured 
by the presence of the factor r(| — s — Ui — v)^, and the same methods as before work to 
give the bound (I3.16p . 
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The same sign case reduces to 

1 \ ^ r r r r r n'^^ m^^ n^^ 

^ \ I I I I / y r /I I w\2t n 



6d,q H^{s,Ui,U2,v,w;Kj)cosh{7rKj) , , , , 
(Xjid) -, — 7 —dvduiduodwds, 



where 



^ w w w 

H_{s, u, V, w; t) = cos(7r(s + Ui + v + — ))r(s + Ui + v + — + it)T[s + Ui + v + — - it) 

£4 

r(i - S - ^2 + f + f ) SW 

Again, we truncate the integrals over s,Ui,U2,'w at imaginary part <^ {qKY . The conver- 
gence in f-aspect can be detected via Stirhng's approximation, which gives 

cos(7r(s + ui + v + Y))r(s + ui + v + — + it)T[s + ui + v + — -it) 



sin(7r(s + Ui + v))V{j2 — s — Ui — v)' 



,V{v)Va-s-U2 + 



T{\-S-U2+V+^) 

A careful but elementary estimation gives the bound 



00 



e "e 



Thus we may take = — ^ — 3£: and Cu^ = h — 3e to get 



2 



'M,N 

Notice this is better than fl3.16p . Actually, this phenomenon holds true for Em+ also: for 
this term, the convergence in kj aspect holds even when all the lines of integration are at e, 
say, and the loss of (N/M)^ does not appear in the bound. 

The treatment of the continuous spectrum follows similar lines. 
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